Superposition in nonlinear wave and evolution 
equations 
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Real and bounded elliptic solutions suitable for applying the Khare-Sukhatme su- 
perposition procedure are presented and used to generate superposition solutions 
of the generalized modified Kadomtsev-Petviashvili equation (gmKPE) and the 
nonlinear cubic-quintic Schrodinger equation (NLCQSE). 
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1 Introduction 



As has been shown recently |Cooper et al. 2002 , |Khare et al. 2002a] . |Khare et al. 2003] . 



Kh are et al. 2002bj (periodic) Jacobian elliptic functions (if they are solutions of 
a certain nonlinear wave and evolution equation (NLWEE)) are start solutions for 
generating new solutions of the NLWEE by a linear superposition procedure. Thus, 
elliptic functions are of specific importance for finding solutions of NLWEEs. On 
the other hand, based on a symmetry reduction method, a technique to obtain ellip- 
tic solutions of certain NLWEEs was proposed and applied to the gmKPE and the 
NLCQSE |Schiirmann 1996], jSchurmann et al. 2004a] . jSchurmann et al. 20 04b . 

department of Physics, University of Osnabriick, D-49069 Osnabriick, Germany, cor- 
responding author: jnickel@uos.de. 

2 Department of Mathematical Sciences, University of Oulu, FIN-90014 Finland. 



1 



It is the aim of the present paper to combine these approaches in order to ob- 
tain general elliptic solutions that can serve as start solutions for superposition 
(" suitable solutions" ) . 

The superposition procedure can be described as follows Cooper et al. 2002] : If a 



solution of a NLWEE[-0(x, y, t)] = can be expressed in terms of Jacobian elliptic 
functions 

I 

^(x, y,t) = a u qn u \p(x + ky + vt), m] , (1) 

where qn is anyone of the Jacobian elliptic functions and a u ,fi,c are constants, 
then the superposition solution Cooper et al. 200*2] Eqs. (4), (14)] 



'" ! ' \ n(X-l)K(m) 

fi(x + ky + Vpt) H , m 



(2) 



A=l u=0 

where n € {2, 4} (depending on the periodicity of the Jacobian elliptic function 
and on v) and K(m),m denote the complete elliptic integral of first kind and 
the modulus parameter (0 < m < 1), respectively, also may be a solution of the 
NLWEE. The number p is integer (it depends on the NLWEE whether it is even or/ 
and odd) and the speed v p can be determined by using certain remarkable, recently 
established, identities involving Jacobian elliptic functions |Khare et al. 2002a] . 
Kh are et al. 2003] . It should be noted, that the existence of solutions ("I} of 
a certain NLWEE does not necessarily imply the existence of a solution (J2J). As 
shown in Refs. |Cooper et al. 2002] , |Khare et al. 2002b] solutions (O exist for the 
Korteweg-de Vries equation (KdV), the Kadomtsev-Petviashvili equation (KP), 
the nonlinear (cubic) Schrodinger equation (NLSE), the Ac/) 4 — field equation, the 
Sine-Gordon equation and the Boussinesq equation. On the other hand, it may 
happen, as will be seen below, that a solution is known but does not lead to 
a solution ©. It is crucial for the procedure, that appropriate relations between 
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Jacobian elliptic functions are known. 



The symmetry reduction approach can be outlined as follows ISchiirmann et al. 2004b] : 
The NLWEE[^(j;, t)] = is reduced by an appropriate transformation tp — ► / 
(e. g., ip(x,t) = f(z),z = x — ct), where / is supposed to obey the ordinary 
nonlinear differential equation ("basic equation") 



2 = af + 4/3 f + 6 7 / 2 + 45/ + e = R(f), (3) 

(with real z, f(z), a, (3, 7, 6, e), leading to an equation P(f) = 0, where P denotes 
a polynomial in /. Vanishing coefficients in the polynomial equation P(f) = im- 
ply equations which partly determine the coefficients a, f3, 7, 5, e in Eq. (0). In 
general, the coefficients depend on the structure and parameters of the NLWEE 
and, finally, on the parameters of the transformation tp — > f. Thus, the problem 
of finding a solution of the NLWEE is reduced to finding an appropriate transfor- 
mation that leads to the basic equation ©. 



As is well known |Weierstrass 1915, Whittaker et al. 1927 the solution f(z) of Eq. 



(J2J) can be written as 

/(*) = /0 + 7fZ7 R '\ f0 \ m(f s V W 

where /o is a simple root of R(f) and the prime denotes differentiation with 
respect to /. 

The invariants 52,53 of Weierstrass' elliptic function p(z;g2,g3) are related to the 
coefficients of R(f) by |Chandrasekharan 1985 

g 2 = ae-4/35 + 3 7 2 , (5) 



3 



g 3 = a~/e + 2/37(5 - a5 2 - 7 3 - e/3 2 . 
The discriminant (of p and R [Chandrasekharan 1985 ) 



(6) 



A = gl - n g \ 



(7) 



is suitable to classify the behaviour of f(z). The conditions 



A ^ or A = 0, 52 > 0, g 3 > 0. 



(8) 



lead to periodic solutions ISc hiirmann et al. 2004b] . whereas the conditions |Abramowitz et al. 1972j 



are associated with solitary wave like solutions. 

Physical solutions f(z) must be real and bounded. Considering the phase dia- 
gram of R(f) |Schiirmann 1996] . |Drazin 1983j one obtains conditions, expressed 
in terms of the coefficients of the basic equation, that determine physical solutions. 
For convenience these conditions are referred to as the phase diagram conditions 
(PDC) in the following. 



2 Elliptic start solutions for superposition 



To apply the superposition procedure it is important to know whether a solution of 
the NLWEE according to (pQ) exists. To check this it is useful to rewrite Weierstrass' 



A = 0, g 2 > 0, g 3 < 



(9) 



function p as [2j 




(10) 



where e\ > e 2 > e% denote the roots of the equation 
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4s 3 _ g2S -g 3 =0. (11) 

Substitution of Eq. (QUI) into Eq. © yields [Sj 



f(z) = (a/o 3 + 4/3/p 2 + 2e 3 / + 5 7 / + 2£)sn 2 (yiT^^, m) + 2(ei - e 3 )/o 
(-«/o 2 - 2/?/ + 2e 3 - 7)sn 2 (^T^^, H + 2(ei - e 3 ) 

(12) 

with m = g^Zeg • Comparison with Eq. Q shows that 



- a/ 2 - 2/3/o + 2e 3 - 7 = (13) 

is a necessary and sufficient condition that defines the subset of solutions (^Q). If 
a = holds the simple root /o of can be choosen such that Eq. (fTT?|) and 

PDC are satisfied [3]. If a 7^ and (3 = 5 = Eq. (|13j) is satisfied also. If 
a j= and /3 / 0, 5 = e = 0, Eq. (JT3J), PDC, and the condition A = 0, g 3 > 
are not consistent, so that trigonometric functions (which are possible for A = 0, 
g 3 > 0) are not suitable for superposition, because f(z) is a constant according to 
the general solution of Eq. (J3J) 1 . 

Equation (|TT?|) represents a relation between the parameters {a, /?, 7, 5, e} and thus 
determines a subset of parameters of the problem modelled by the NLWEE for 
which further solutions can be generated by superposition according to Eq. |2J). 
Combining Eqs. (|12|) and (|13|) (with a = 0) we obtain 

£ , \ 2e 3 - 7 12e§ - 3 7 2 + 4(35 2 . , . 

/(*) = + 4 3 ^ (ei _ e3) ^(ViT^^, m) , (14) 

where ei, e 3 must be chosen as the largest and smallest roots of Eq. (|11|). re- 
spectively, so that the condition (fT3)) is valid for a simple root /o of Eq. © that 
satisfies the PDC. 
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Equation (fTl)) can be evaluated explicitely subject to the two cases a = and 
a ^ 0, (3 = 5 = 0, respectively. If a = and, for simplicity, e = the start 
solutions for superposition are 



15 



-37+A/97 2 -16/3<5 k 



sn 2 W-3 7 + \/ 9 7 2 - 16/3<5*, 



37+y / 97 2 - 16/3.5 
37- a/97 2 -16/35 



, /ft > 0, 7 < 0, 



37+ a/97 2 -16/35 ^ 2 ( (9 7 2 - 16/35)1 37+^972-16/35 

Z ' 2^972-16/35 



4/3 



-cn 



06 < 0, 
(15) 

where the various possibilities to satisfy (|11[) and (|13j) have been taken into ac- 
count and A = 4/3 2 <5 2 (97 2 — 16/35) > is necessary and sufficient to fulfill PDC 
[Bronstein eTa^ OOO . 



If a 7^ 0, (3 = 5 = the start solutions read 



37+ 



Hz) 



|^dn 2 J37 + y/^^Tez, 



2^ 



37+v / 97 2: 



-37+^/97^ 



sn 2 I y— 37 + -y/97 2 — aez, 



37+\/9t 2: 



37- 



, a < 0, 7 > 0, e < 0, 



, a > 0, 7 < 0, e > 0, 



_ 37+ V97 2 -ae cn2 / ^g,^ _ .1^ 3tW5£=^\ 



a < 0, e > 0, 
(16) 

where A = 64a 2 e 2 (97 2 — ae) > and - according to the Cartesian sign rule - three 
numbers of sign reversals in the sequence of coefficients of R(h) or A > and 
a > and two sign reversals to fulfill PDC. 
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To sum up, Eqs. (fT5|) and <|16|) represent all elliptic solutions with a = 0, e = and 
a ^ 0, f3 = 6 = 0, respectively, that are suitable for the the procedure suggested 
by Khare and Sukhatme. "All elliptic" means that the solutions presented in Refs. 
Cooper et al. 2002| , [Kb" are et al. 2002b| are particular cases of Eqs. ()15[). (|16[). 



"Suitable" includes that the superposition procedure may fail if solutions according 
to Eq. (|15|) or (J16|) are inserted into the NLWEE in question leading to conditions 
that cannot be evaluated with respect to v p (cf. Eq. ©) because the associated 
relations between Jacobian functions are unknown (cf. Sec. EJ). Examples to 
obtain superposition solutions are presented in the following. Equation ()14j) can 
be evaluated in the same manner subject to the PDC to obtain physical elliptic 
solutions if the simplifying assumption e = does not hold. 



3 Superposition solutions of the generalized 
modified Kadomtsev-Petviashvili equation 

The approach outlined in the previous section can be elucidated by investigation 
of the gmKPE (Superposition solutions of the NLCQSE are presented in EI) 

Tpxt + ((a + b^) ^ x ) x + c^ xxxx - o 2 4> yy = 0, (17) 

where a, b, c, q, a 2 are real constants. As shown previously [Schiirmann et al. 2004a] 
elliptic traveling- wave solutions to Eq. ()17|) exist. The set of these solutions is de- 
termined by 



ip(x,y,t) 



z 

(■2 



x + ky + vt, 

af + 4(3 f + 6 7 / 2 + 45 f + e, 



(18) 



where a, /?, 7, 5, e are given by Eqs. (16a)-(16g) in Ref. |Schiirmann et al. 2004a] . 
As shown above, the conditions a = or (3 = 5 = 0, a ^ lead to suitable 
solutions. Imposing additionally the PDC and condition respectively, the 

parameters of solutions (|18|) are 



q = \,a = -^- c ,P = 0n= k ^ V - ,S = 0,e^0,c^0, (19) 

a, /c 2 2 

g = l,« = 0,/3 = -^, 7 = a Q ~ V ,e = 0,c^0, (20) 

r t 2 2 

q= l,a = = 0,7 = — ,«5 = 0,e / 0,c ^ 0, (21) 

DC DC 

g = 2 ,a = 0,/3 = -^, 7 = 2(fc2 ^~ 7j) , e = 0,c^0.. (22) 

Referring to (|T§)) and (|2*U)) first, solutions according to (fTB|) and (|15f) . respectively, 
have to be evaluated. Inserting ()19|) into (|16|) . one obtains the suitable start 
solutions 



ij){x,y,t) = < 



Bj dn^ [ W (x + ky + ut), mi] , | > 0, > 0, e < 



B 2 sn 2 [/i 2 (x + ky + vt), m 2 ] , £ < 0, < 0, e > 0, (23) 



B 3 cn 2 [/i 3 (x + ky + vt), m 3 ] , \ > 0, e > 0, 

where Bj, fij, rrij are determined by inserting the parameters a, 7, e according to 
()19l) into (|16|) . Formally the same result is obtained by inserting (|20|) into (|15l) . 
Referring, secondly, to ()21() the solutions follow from (|16|) as 
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B 1 dn + ky + irt),mi] , | > 0, > 0, e < 0, 



B 2 sn [ M2 (x + ky + vt),m 2 ] J < 0, < 0, e > 0, (24) 



5 3 cn [/i 3 (x + ky + vt), m 3 ] , | > 0, e > 0, 

where (again) Uj, /i,-, mj are determined from (|16|) with parameters according to 
H21j) . Formally the same results are given by (|22|) and (|15|) . 

According to Eq. © the first solution in ()24|) leads to a superposition solution for 
p = 2 



ip(x,y,t) = B dn (//(x + fcy + t> 2 t) + (i — l)K(m), m) 



i=i 



B = Bi,fi = Hi, m = mi. 



(25) 



Inserting ^(x, y, t) (denoting dj = dn (fJ,(x + ky + v 2 t) + (i — l)K(m),m)) into Eq. 
((T7|) (a = 0, because f3 = according to l|21j)l we get 



d 2 d 2 

(—Bmuv 2 — Bc/j 3 (2m — m 2 )) — SjCj + o 2 Bmuk— SjCj 

Ht. — i dy ^ — ' 



8=1 



i=l 



d 



+2bB 3 m 2 fj? Yj di [Yj SiCi j - m ^ bB3 [12 di J J2 SiCi ( 26 ) 



1=1 \i=l 



d 



+6Bcmu 3 — d. 2 SjCj = 0. 
dx ' 



i=l 



The last three terms of Eq. (|2fi[) can be simplified as follows. 

Using did 2 = y/l — m and cisid 2 + c 2 s 2 di = |Khare et nl 2n02a1 Eqs. (31), 
(39)] we obtain 
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X^ di SiCi = ^i s i c i + Vl - m ^ SjCj. (27) 

ten as 



, 2 

{-Bmu,V2 - Bcu 3 (2m - m 2 ) - mubB^yJl - m) — SjCj (28) 
v dx ^-^ 

i=l 

d 2 d 2 

+a 2 Bm^k— \^ SjCj + 8m/j (6c/U 2 — 6-B 2 ) — dfsjCj = 0. 

y i=l i=i 
The expression (6cfj?—bB 2 ) vanishes identically jH). With ^ ]Ci=i s « c i = ^ S 2 =i s i c « 
Eq. |(2Hl reads 



d 2 

(—Bmnv2 — Bcu?(2m — m 2 ) — mubB 3 -\/l — m + a 2 Bmuk 2 ) — — > SjCj = 0, 

i=i 

(29) 

so that the speed V2 is given by 



V2 = -c/x 2 (2 - m) - bB 2 y/l - m + a 2 k 2 . (30) 

Thus, we have found a superposition solution of Eq. 1|17|) for this particular case. 
The start solution and the superposition solution are shown in Fig. ^ 
We can generate superposition solutions for p = 3 from (J23|). As an example we 
consider the solution of the form dn 2 in detail and compare it with the results 
of Cooper, Khare and Sukhatme Cooper et al. 2002| . According to Eq. (J2J) the 



superposition ansatz reads 



$(x,y,t) = B ^2 dn2 (v( x + k V + v 3t) + — ; m ^ 
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Figure 1: The start solution ip(z) (cf. first solution of Eq. (|24|)) and the super- 
position solution tp(z) (cf. Eq. (|25|0 for a = —2, 7 = 4, e = — 1, z = x + ky + vt 
and 2 = x + % + v^t, respectively. 



.B = B%, /j, = hi, m = mi. 



(31) 



Inserting tp(x, y, t) (denoting dj = dn ( fi(x + ky + t>3t) + 2 ^ 1 )^"( m ) ^ m j ) j n ^ jjq_ 
((T7I) . we obtain 



d 3 ^ 3 

2Bmu(v3 + 8cu 2 — Acmu 2 ) — djSjCj + 2a 2 Bmiik— djSjCj 



i=i y i=i 

\ 2 3,3 

a 



-AB 2 bm 2 fi 2 d i s i c i ~ ZB 2 bmn ^ d 2 — ^ d^ (32) 

\j = l / j=l X 8=1 



+24 J Bcm// 3 — y^dfsjCj = 0. 



i=i 



The last three terms can be rewritten as 



' 3 



3 , 3 

a 



- 2B 2 bmfi \ -2mfi [ ^djSjCj ) + ^ d 2 ^ ^ djSjCi - 12 ^"^ d ? SiCi ) ' 



vi=l 



i=l i=l 



i=l 



(33) 
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whereas evaluation of fSi=i Ylj^i ^j s j c i) yields 



3 \ / 3 \ ^ 3 3 

\ Z d * Yl d « = - 2m / i ( Z diSiCi ) + Z d2 ^ Y d ^ ci 

3 

.. , , * — * / , (l \ 

(34) 



+2m^^d 2 s 2 c 2 - U^diSiC-) 
i=i i=i ^ 1 ' 



( 3 \ 2 3 d 3 <i 3 

-2m/i ( ^ diSidj + Z ^ ^ Z d iSjCi - — ^ df SiCi 



2 

Because 12|g- = 1 (in Eq. (JHSJ) holds identically, we can use Eq. (|34j) and 
[Khare et a/. 2002 bl Eq. (11)] 



d I 3 \ d 3 

^ ( Z d ? Z d i s i c i = ^( 3 > m )fa Y diSiCi ' ( 35 ) 

i=l j^i I i=l 



to rewrite Eq. (|3*2*|) as 



d , „ o „ , d 



-2Smu(w3 + 8cu 2 -4cmu 2 + £?&yl(3, m)) — djSjCj+2cr 2 Bnifik — djSjCj = 0. 

dx U d yU 

(36) 



2Bm/i(u 3 + 8c^ 2 - 4cm<u 2 - a 2 k 2 + BbA(3, m))-^- V d^c; = 0. (37) 



Using ^ i=1 djSjCj = A;^ £)i=i djSjCj this equation reads 

3 

(fa 

Thus, the speed i>3 in the superposition solution (|3~Tj) (of a particular case) of Eq. 
(fT7|) is given by 

v 3 = 4cm/x 2 + a 2 k 2 - 8c/j 2 - BbA{3, m). (38) 
The start solution and the superposition solution are shown in Fig. [21 



12 




1 2 3 4 5 6 7 



Figure 2: The start solution ip(z) (cf. first solution of Eq. I)23[l) and the super- 
position solution ip(z) (cf. Eq. (|31jll for a = —1, 7 = 1, e = —1, z = x + % + u t 
and z = rr + ky + i^t, respectively. 



Applying an analogous procedure with the ansatz 

7, \ 2 ( / , v 2(»- l)if(m) \ 
^(x,y,t) = B 2^ sn 2 (^(x + ky + w 3 t) + -i g ; , m j , 

B = B 2 ,n = H2,m = m 2 (39) 
and with the ansatz 

A 2 / 2 (2 l)if(m) \ 

v(x,y,t) = B 2^ cn ^(2; + «y + ^3*) H g 

■B = B 3 ,fi = fi 3 ,m = m 3 (40) 
we obtain superposition solutions with 

^3 = 4cm M 2 + 4c/i 2 + a 2 A; 2 + fl& A(3 ' m) ~ 2 (41) 

m 

for solution (|31?)) and 
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v 3 = -ScrV + 4c/x 2 + a 2 k 2 - Bb A ^ m) 2(1 m) (42) 

m 

for solution (|4"0|). 

In deriving 1)4 1JI and l)42j) we have used the relations 

d f 3 3 \ 1 d 3 

^ E s ? E s i c i d i = w> m ) - 2 ) Tx E s » c » d * ( 43 ) 

\*=1 } i=1 



and 



d ( 3 3 \ 1 <i 3 

\i=i j^i J i=i 

respectively, which follow from Eq. ()35jl and well known relations between Jaco- 
bian elliptic functions. 



Comparison of the Kadomtsev-Petivashvili equation together with the ansatz con- 
sidered by Cooper, Khare and Sukhatme Cooper et al. 2002] Eqs. (l),(4)]with the 
Eqs. (|T7j) . (fT§)) and (f3"H) shows that, apart from an additive constant |Jaworski et al. 2003| . 
our result ()38|) is consistent with that of Cooper, Khare and Sukhatme |Cooper et al. 2002 



Eq. (11), P = 0] .The cases related to (|4"Tj) . (|l2*)l have not been considered in Ref. 
Cooper et al. 2002 . 



To conclude, we note that real and bounded suitable solutions of the gmKPE only 
exist for four different values of q (cf. 1)191) - (|22|0 . though there is no restric- 
tion for q (apart from being real) of the known elliptic solutions of the gmKPE 
Sch urmann et al. 2004aj . 

The second of Eqs. does not lead to a superposition solution although the 

solution has the form (^Q) 6 . In this case, it seems that an appropriate identity 
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for Jacobian elliptic functions does not exist. Thus, the claim at the end of Ref. 
|Cooper et al. 2002] seems to strong. 



4 Summary and concluding remarks 



By combining the superposition principle and symmetry reduction we obtained 
general elliptic solutions suitable for superposition. The results were applied to 
the gmKPE and the NLCQSE (seeEJ). In Ref. |Cooper et al. 2002| particular 
elliptic solutions for generating superposition solutions of the NLSE and the KPE 
were used. As outlined above we start from (general) suitable solutions (cf. Eqs. 
(|15jl. Qlrij). (|2Hjl. (J2U)) to obtain superposition solutions more general than those 
in Ref. |Cooper et al. 2002 . We note that there are no restrictions in advance for 
the coefficients of the NLSE and the KPE. Constraints result from the condition 
that suitable solutions exist (cf. Eq. and from the PDC. As is obvious from 

the following list [I] there are rather many NLWEEs that exhibit suitable elliptic 
solutions. Thus, it seems interesting to check whether they lead to superposition 
solutions by applying Eqs. (|15j) and (fTrT)) . 



Table I: Elliptic solutions of various nonlinear wave and evo- 
lution equations. 



Equation 


Ansatz 


Basic equation 


Suitable for 
superposition 


tpt ~ 1p1px - Ipxxt = 

Benjamin-Bona- 
Mahony 


V> = f(kx - ct) = f(z) 


(f\ 2 - f Z + f' 2 
+45 f + e 


+ 


i>tt ~ ipxx + 3(ip 2 ) xx 

Ipxxxx — 


V> = f{kx - ct) = f(z) 


( t \2 _ 2/ 3 . c 2 -k 2 f 2 
\Jz) — -%r T fc 4 J 

+45 f + e 


+ 
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Boussinesq 








i>t + ipi>x ~ ipxx = o 


V> = f{kx - ct) = f(z) 


/ f \2 (,^c / — fcc/ +4fc 




Burgers 








tptt ~ ipxx + sin ip 


-0 = 4 arctan[/(/ca; — ct)] 


(/ 2 ) 2 =3( 7 + I (^ T )/ 4 




+isin| = 


= 4 arctan[/(z)] 


+67/ 2 + 3 7 +(^, 


+ 


Double sine-Gordon 




c 2 / A; 2 




tpt + a^xx ~ H 


v> = 


WsJ 3dai+a 2 (2-d 2 )-^ 




-c | V | 2 ^ = 


f(x)exp(ig(x))exp(i\t), 






a = ai+ia 2 , b = bi+ib 2 , 


g x (x) = dff 


/ A-|- — 62 \ j?2 ISl 
+ (^a 2 -a 2 -2daj / l^J 


+ 


c = ci + ic 2 








Ginzburg-Landau 








^t+{H+^x+^ttx = o 


ip(x,t) = 


( f \2 _ f 3 1 c-fc j-2 




Joseph-hgri 


j(kx - ct) = i(z) 


2Ci j- _ 2C 2 
c 2 k c 2 k ' 


+ 


jHereman et al. 1986J 




Ci, C2 const. 




ipt + a^ 2 ^ + bip x tp xx + 


ip(x,t) = 


(/ z ) 2 = 4/3/ 3 + 6 7 / 2 




S^P^ftxxx ~i~ Ipxxxx — 


f(x - ct) = f{z) 


+45/ + e 


+ 


Korteweg-de Vries 








4>t - Ipxx - Glplpx + Ipxxx 


i; = f(kx-ct) = f 2 (z) 


2 / 2 




= 




/ f 2 of 1 \ 

. 1 1— 4- 2]- + J_ ) 




Korteweg-de Vries- 








Burgers 








ipt + a^x + ipxxx 


%p{x,y,z,t) = /(£) 


(A) 2 = -sf>/ 3 




~\~'4 , xyy ~\~ Ipxzz — 


£ = kx + ly + mz + ujt 


_ ^ f 2 _ 2Ci 


+ 


KdV-Zakharov- 




p 2 = /c 2 + / 2 + m 2 , 




Kusnetzov 








[Baldwin et al 2004 




Ci const. 
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V t = d 3 V + d 3 V 
+3d(uV) + 3d(uV), 
du = dV 
Novikov-Veselov 


V(x,y,t)=ip(z) 
z = x + ky — vt 


( f \2 _ 8o 2 f3 
\Jz) l+k' 2 J 

24a f2 | 12F f2 
l+k 2 -> 1 3(3fc 2 -l)-' 

+W + e, 

F = v + 3C + 3kd; 
Co, C\ const. 


+ 


ipxx - iptt - sin ip = 
sine- Gordon 


ifj(x, t) = 4 arctan[ : ^y] 


= aX A + 67X 2 + e 

(fi) a = -RiCn 

= aT 4 + (67 - 1)T 2 - e 


+ 


-fJ-ipxt = 
Wadati, Segur, 
Ablowitz 

[Wadati et at. 1992| 


ip(x,t) = f{z)e i{ - rx - x ^ 
z = kx — ct 


C f l 2 — CT f 4 

c(c+kfi) J 
c 2 /i(c+fc/i) J 

A' const - 


+ 



If condition 1)13(1 can be fulfilled (e. g., by choosing an appro- 
priate constant of integration) start solutions for superposi- 
tion can be obtained by Eqs. (|15|). 1)16(1 (marked by a "+", 
otherwise by a "-"). 
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A Superposition solutions of the nonlinear 

cubic-quintic Schrodinger equation (NLCQSE) 

Following the lines described in Sees. ^|2]the NLCQSE 

^ + fe-(9i|Vf + 92|V'| 4 )V' = 0, (45) 
(<Zii Q2 real constants) can be solved by applying the transformation 

ijj(x, t) = f(z) exp[i(At + r(z))], z = x - ct. (46) 
Separating real and imaginary parts, we obtain 

qif(z) 3 + Q2f(z) 5 - f"(z) + f(z)(X - cr'(z) + r'(zf) = 0, (47) 

f'(z)(c-2r'(z))-f(z)r"(z)=Q, (48) 

where the prime denotes differentiation with respect to z. 
Equation (|48j) can be integrated to yield 

with C\ a constant of integration. 

Inserting r'{z) into Eq. 1)471) and integrating the resulting expression leads to an 
equation where h = f 2 can be introduced. Thus, we find a basic equation R(h) 
(cf. Eq. (j3J), / — > h) with the following coefficients: 

4 1 4A — r 2 

a = -q2, (3=-qi,l=——,S = 2C 2 ,e = -4C 2 1 , (50) 
6 2 b 
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where C2 is a constant of integration. 

If Q2 = and C\ = all physical solutions suitable for superposition are represented 
by Eqs. (|15j) (/ — > /i). The superposition solutions for p = 3 are given by (cf. Eqs. 



,1 



4(t-l)A'(m) 



exp{i [Xt + (x - v 3 t)f]} 



4(A - /j 2 (2ml(3, m) + (2m - 1))), 



(51) 



ip(x,t) =a^ =1 dn (i(x-v 3 t) + 



4 



2(i-l)K(m) 



m 



exp {i [Xt + (x — v 3 t)^f] } , 



4(A + fi 2 (m - 2) - aW(3, m)), 



(52) 



*l>{x,t) =al]^ = isn ^(x - i> 3 t) + 4( * jHAH ; m exp |i [At + (x - v 3 t)^] } , 
u| =4(A + / u 2 (m + l) + 2ma^V(3,m)). (53) 

To evaluate the speed v 3 we have used in Ref. |Cooper et al. 2002 the Eqs. (8), 

(70), (72), Eqs. (8), (66), (68) and Eqs. (8), (57), (59), respectively. 

It should be mentioned that the start solutions (fH))) suitable for superposition are 



consistent with those of Cooper, Khare and Sukhatme Cooper et al. 2002 . Nev- 
ertheless, the speed ^3 according to Eqs. (|STj) . (|52|). (f53*|) is not identical with v 3 
according to Eqs. (33), (28), (45) in Ref. |Cooper et al. 2002| . Thus, the superpo- 
sition solutions are not determined uniquely. Different identities between Jacobian 
elliptic functions used lead to (in general) different superposition solutions. Ap- 
plying the procedure outlined in Sec. El if c/2 7^ (a ^ 0), f3 = 5 = 0, e arbitrary, 
PDC implies either 52 = (a = 0) or C\ = (e = 0). The choice q2 = (in 
addition to q± = ((3 = 0)) is not of interest, because it leads to a linear equation 
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(|45|). For Cf = we obtain solitary traveling- waves. Thus, since if>(x,t) is not 
periodic, superposition solutions are not possible in this case. 
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[1] The general solution of Eq. (j2J) reads 

IWeierstrass 1915, Whittaker et al. 1927 

f()=f yj^Mf|M3) + i#(/ )[p(, ;g2 , g3 ) - ±R"(f )] + ±R(f )R'"(f ) 
[Z) 2[p(z;g 2 ,g 3 ) - ±R"(f )] 2 - ^R(fo)R""(fo) 
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[2] We assumed A > 0. If A < 0, substitution of p(z) according to Ref. 
[Abramowitz et al. 1972j . 18.9.11 does not lead to an expression of form (^Q). 
Furthermore, as will be seen below, we obtain a polynomial R(f) or R(h), 
h = f 2 , of third degree for generating new solutions by linear superposition. 
But the PDC is not satisfied for a third-degree polynomial with A < 0, be- 
cause two of the three roots are complex conjugate. 

[3] Since we are interested in physical periodic solutions we can always assume 
that a simple root exists. 

[4] If (3 = holds 7 must be negative otherwise Eq. (|13|) and PDC are not fulfilled. 
For /o to be a simple root 5 2 — |e7 must be positive. If (5 ^ the discriminant 
A does not vanish, so that /o = e |J 7 is a simple root of R(f). 

[5] If parameters according to Eqs. (|21|) are inserted into Eq. Q16JI one obtains B, 
fjL, m so that 6-Bcm/i 3 — m^bB s in Eq. (|28j) vanishes identically. 

[6] An ansatz ifj = cn or ip = sn (cf. eq. 1)25(1 ) leads to equations which have the 
form (J26|) . Because there is no relation C1C2 = const, and S1S2 = const., respec- 
tively, there is no possibility to replace the appearing sums ^ Cj djSj) 2 , 

E c «) 2 5E d « s « and E s i (E c i d i) 2 > (J2 s if 7E E c idi, respectively. Up to 
our knowledge there is no appropriate relation involving Jacobian elliptic 
functions that would simplify the equations similar to (|2*^|) . so that the speed 
V2 for which t/j is a solution of eq. (|17|) cannot be determined. 

[7] As outlined above start solutions for linear superposition can be obtained if 
a = Oor/3 = <5 = 0. In the case of the Burgers equation these conditions lead 
to k — > oo or c = for a traveling wave ansatz ip = f(kx — ct). 

[8] d only depends on the coefficients of the Ginzburg-Landau equation a, b, c. 
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[9] As outlined above start solutions for linear superposition can be obtained if 
a = or (3 = 5 = 0. In the KdV-Burgers equation the only parameter which 
can be varied in the basic equation is k. If k — > oo then a — > 0, but also 
(/ z ) 2 — > so that / = const.. 
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